UNCLASSIFIED 


AD NUMBER 


LIMITATION CHANGES 
TO: 


Approved for public release; distribution is 
unlimited. Document partially illegible. 


FROM: 


Distribution authorized to U.S. Gov't. agencies 
and their contractors; Critical Technology; MAY 
1965. Other requests shall be referred to Air 
Force Technical Applications Center, 
Washington, DC. Document partially illegible. 


This document contains export-controlled 
technical data. 


AUTHORITY 
usaf ltr 25 jan 1972 


THIS PAGE IS UNCLASSIFIED 


340) BHiILOH ROAD 


PE eT Ie PY HES OTE 


THE BEOQOTET 


INTERPOLATED JEFFREYS AND PULLEN SEISMOLOGICAL TABLES 


ee oe 
ee a eee tr) 


ake a 
Le apr 20 1956 \ 
poe TT 


THIS DOCUMENT 1S SUBIECT TO SPECIAL 
EA ual CONTROLS AND EACH TRANSNITEAL 
TO FOREIGN GOVERNMENTS CR FOREIGN 

WAT iONA L MAY BE MADE GNLY WITH PRIGHK 


APPROVAL OF CHIEF,APIAG, 


HNICAL CORPORATION 


GARLAND, TEXAD 


acvocevssuecsee-el se Ses}. shpheqoG”-aei>ERVaPIOL I DENENOOITUIEASOOFIDOISINNSCSPRONETIIETRSREONOGOOON, 3 RSPOSIRNRETOUSTOONETSTT 


BEST 
AVAILABLE COPY 


TECHNICAL REPORT NO, 65-35 


. INTERPOLATED JEFFREYS AND BULLEN 
SEISMOLOGICAL TABLES 


‘ey 


H. S. Travis 


etn SSA SHC ayTET HC IpeCTTSERIeboncteSeromumetearwetionin Fy eegrnier sererrr 


AFTAC Project No: VT /40° * 
ARPA Oxder No: 104.-61 
ARPA Code No: 8100 
. Contractor: The G.otechnical Corporation, 
Garland, Texas 
Contract No: AF 33(657)-12145 


coor be aaa ON 


THE GECTECHNICAL CORPORATION 
3401 Skhilch Road TS PICU INT 16 SUBIECT TG SPECIAL 
2 Garland, Texas , is ard ZACH TRANEETTAL 
15e NMZAYS OR FOREN 
WATIQUAL MAY BE MADE ONLY WITH PRIOR 
APPHOVAL OF CHIEF AFTAG. 
5 May 1965 
i 


SONTENTS: 


Page 

ABSTRACT 

i INTRODUCTION 1 | 
i 

2. DESCRIPTION OF METHOD 2 | 
3, NUMERICAL TECHNIQU: x 
3.1 Lagrangian interpolation formula 3 i 

3.2 Linear interpolation formula 4 

4, REMARKS 4 | 
5. INDEX OF PHASES 5 i 
6. INTERPOLATED JEFFREYS AND BULLEN 
SEISMOLOGICAL TABLES 

i 


TR 65-35 


ABSTRACT 


The travel times for all body phases tabulated in the Jeffreys and Bullen 
Seismo]ogical Tables (1958) have been calculated by interpolation for distance 
intervais of 1 degree at each of the 14 focal depths established by Jeffreys and 
Bullen. Travel times have been similarly derived for the phases "diffracted P" 
(P at distances beyond 105 degrees) and "PKPPKPPKP," based on data extrac- 
ted from the Gutenberg and Richter travel-time charts (1934-1939). The cal- 
culations were effected by means of a digital computer program principally 
employing Lagrange's four-point interpolation formula. This report discusses 
the techniques invoived and presents the fully interpolated tables in a convenient 
form. 
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INTERPOLATED JEFFREYS AND BULLEN 
SEISMOLOGICAL TABLES 


1, INTRODUCTION 


This report presents a fully interpolated version of the Jeffreys and Bullen 

‘J-B} Seismological Tables” and discusses the manner in which it was gen- 
erated. The travel times for all teleseismic body phases contained in the 

J-B tables are given for distance intervals of 1 degree at each of the 14 
standard J-B focal depths. The times are tabulated in a convenient reference 
form very similar to that empicyed for tne "P" and "S'' phases in the existing 
J-Btables. Travel times are similarly presented for the phases "diffracted P" 
(P at distances beyond 105 degrees) and "PKPPKPPEF , " based on data extracted 
from the Gutenberg and Richter travel-time charts. 


These interpolated I-B tables have been calculated as a corvenience for seis- 
mologists who must have frequent access to accurate, closely resolved (in 
distance and depth) travel times for all phases. The only phases in the exist- 
ing J-B tables listed for distance intervals of 1 degree at all depths are ''P"' 
and 'S." lt is anticipated that the tables will save m=ny hours of needless, 
repetitious manual interpolation on the part of analysts trying to identify and 
associate phases. Particular application is seen for compiling seismological 
bulletins, sucn as those assembled under the Long-Range Seismic Measure- 
ments Program, which routinely involve a vast amount of table lookup for 
accurate association of phases. 


Section 2 of this report discusses in detail the interpolation and computational 
methods used to create the interpolated tables, and section 3 specifies 
mathematically the interpolation formulas used. Section 5 provides an index 
of the phase3 available, their various branches, and the ranges of distance 


lteffreys, Harold, and Bullen, K. E., 1958, Seismological tables: British 
Association for the Advancement of Science Publication 


“Gutenberg and Richter, 1934-1939, Seismic wave papers, I-IV: Gerlands 
Beitrage Zur Geophysik 
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over which they are tabulated. Depth phases suck as "'pP'' and "'sS"' are not 
included because of the ease with which their times can be calculated from 
the 'P'' and "'S'' tables. 


2. DESCRIPTION OF METHOD 


A program for the Control Data 160-A digita) computer was written to inter- 
polate the existing Jeffreys and Bullen seisrr:ological travel-time tables and 
tabulate the results. Primarily, the program employed the Lagrangian four- 
point interpolation polynomial to interpoiate between the times; however, . 

when four points csuid not be obtained, a linear interpolation method was used. 
For each phase where it was necessary, the program first calculated by inter- 
polation the depth allowances at all standard J-B depths for distance intervals 
of 1 degree. Depth allowances in the Jeffreys and Bullen tables are generally 
given at distance intervals of only 10 degrees. Subsequently, the interpolated 
depth allowances were subtracted from the surface times of the phase to obtain 
the corrected travel times at depth. In some cases, the times for surface focus 
were given at 5-degree intervals and also had to be interpolated for each degree. 


For core pnases other than "PKP," "SKS," "PcP," and ''ScS,'' the allowances 
for depth had to be obtained by first determining the dt/dA at 1-degree intervals 
for each phase's surface travel-time curve. Then the allowanccs fur depth 
were taken as equal to those given (or interpolated) for the above core phases 

at distances where dt/dA were the same. For example, to obtain the. desired 
depth allowance for the phase "SKKS"' at 150 degrees for depth 9. 06R, it is 
determined that dt/dA at that distance is 5.5 sec/degree, which is matched by 
"SKS" at 94 degrees. Therefore, the depth allowance for "SKKS" at 150 degrees 
for depth 0.06R is equal to that for 'SKS"' at 94 degrees and depth 0. 06R. 3 
those phases whose depth allowances were not defined by the J-B tables, such 


as 'PP'' at 10 degrees and depth 0.06R, the times were lefi blank. 


Travel-time tables were tabulated for those phases which produce multiple com- 
binations, such as 'PcSPKP.'' Tables for the phases 'ScPPKP" and 'PKPScP" 
(or PKPPcS} were tabulated by obtaining the prover depth allowances (for ScP 
and PKP, respectively) by the standard procedure and subtracting these times 


3 Jeffreys and Bullen, op. cit. 
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from the surface times of 'PcSPKP." The travel times fer diffracted 'P" 


and 'PKPPKPPKP" were obtained from the Gutenberg and Richter travei- 
time charts at distance intervals of approximately 0.0 degree and at depth 
intervals of 100 kilometers. From these, travel times were determined for 
distance intervals of ] degree at the standard J-5B focal depths. 


Generally, the times for each phase given in the standard Jeffreys and Bullen 
travel-time tables represent average times with standard errors of under 

1 second. # Inevitably, the interpclated tables presented ..erein have the same 
error; however, it should be po!nted out that possibilities of smal! disagree-~- 
ments with the standard tables arise from the smoothing technique introduced. 
In the majority of cases, the interpolated times agree to 0.1 second with the 
standard tables. However, the derivatives (dt/da) of the travel-time curves 
for several of the more exotic body phases (e.g., PKPPKS} were so erratic 
the interpolation technique employed smoothed over these points. Conse- 
quently, the times differ occasionally from a linearly interpolated time. The 
magnitude of the disagreement is generally less than half a second. 


3. NUMERICAL TECHNIQUES 


3.] LAGRANGIAN INTERPOLATION FORMULA 


Numerical four-point interpolation (Hildebrand, F. B., 1956) can be obtained 
from the Lagrangian interpolation polynomia of degree 3 expressed in the 
following general form: 


i) 


u(x) = by 4 Ax) fix) 
k=9 


where: 


y(x) is the Lagrangian interpolation polynomial of degree 3 tnai 
approximates the assumed continuous function f(x). 


£(x},) is a series of discrete ordinate points (corresponding to 


abcissa values xg, X1, X2) and x3) representing the function to 
be approximated. 
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and the desired Lagrangian coefficients are given by: 


(x-xQ) .-- (x-%, 01) (x-xE 1)... (X-%3) 
tax) = = 
(Xp -x Qe +» (KH. TM KKR yy). - » (X_-H3) 


A eS Ra a ere ame 


where che factors (x-x,) and {x;,-x:) are to be omitted in the products. Thus, 
avaJv.. of y may be calculated for any value of x by means of the Lagrangian 
polynomial. To obtain the highest degree of accuracy, the value of x should 
be centered between the two innermost data pcinte (x, and x2). 
3.2 LINEAR INTERPOLATION FORMULA 
Linear numerical interpolation (Hildebrand, F. B., 1956) between two tabular 
points can be obtained from the following formula: 

fx) = flxg) + PFO | aise) - x9) 

(x1 -xQ9) 

where: 


f(x) is a function assumed to be linear between the two 
tabular points, xg and xj,. 


4, REMARKS 


The intervolated Jeffreys and Bullen seismological tables are also available 
on punched cards. 
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PHASE 

l. P 

2. P diff 

3S 

4, PP 

52 PS 

6. SS 

%. SP 

8. PPP 

9. PPS 
10. PSS 
11. SSS 
12. SSP 
13. SPP 
14. FecP 
15. PcS 
16. ScS 
lve ScP 
18. PKP 
i9. PKS 
20. SKS 
21. SKP 
22. PKKP 
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5. INDEX OF PHASES 


BRANCH 


AB 
BC 
DEF 


AB 
BG 
DEF 


AC 
DEF 


AB 
BC 
DEF 


AB 
BC 
DF 


143° -180° 
143° -147° 
110° -180° 


13 )° - 148° 
139° ~140° 
Poe? -18e 


70-133" 
99° -180° 


130° -148° 
13G° -149° 
14° -180° 


104° -126° 
93° ~126° 
0” = 162" 


OT ae a 


76 


79 


PHASE 


23. 


24, 


z5. 


26. 


27. 


30. 


cg 


32. 


ce 


34. 


Er, 


36. 


PKKS 


SKKS 


SKKKS 


PoPRiP 

1 acu 
ScoP 
ScSPKP 
ScPPKP 
BKPP or 
PKFPcS or 
PRIPScP 

Ss aa 
PREPRESS 


PKPPKS or 
PKPSKP 
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DF 
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DF 
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DISTANCE xANGES 


134° - 146° 
100° -146 
0°-168° 


85° -243° 
o” -173° 


134° -146° 
100° ~- 146° 

0° -168° 
90° -370° 


82° -187° 
132° -186" 


126° - 194° 
127°-180° 


112°-140° 
112° -230° 


160° -206° 
121°-180° 


110°-194° 
127 °-180° 


82° -187" 
132° -180 


110°-194° 
127° -162° 


160°-206° 
121°-180° 


0°-140° 


75° = ge" 
0°-145° 
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CORES HAIER 


PHASE BRANCH DISTANCE RANGES PAGE 
37. PKSPKP 
B 15 - 84° 150 
wDF 0° ~145° 150 
38. SKPPKP 
B dna re 154 
DF 6G =145 154 
39. SKSP 115° -205° 158 
40. PKPPKI-PKP 44° ~104° 161 
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